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We present a semiclassical theory for electron drag between two parallel 
two-dimensional electron systems in a strong magnetic field, which provides a 
transparent picture of the most salient qualitative features of anomalous drag 
phenomena observed in recent experiments, especially the striking sign reversal of 
drag at mismatched densities. The sign of the drag is determined by the curvature 
of the effective dispersion relation obeyed by the drift motion of the electrons in 
a smooth disorder potential. Localization plays a role in explaining activated low 
temperature behavior, but is not crucial for anomalous drag per se. 
PACS: 73.21.-b, 73.63.-b, 73.43.-f 



Electron drag in double-layer two-dimensional electron systems has been established as 
a valuable probe of the electronic state within each layer, and also of interlayer interactions 
Q]. In a drag experiment a current is driven through one of the layers, which, via interlayer 
scattering, produces a drag voltage in the other layer. Usually the drag is positive in 
the sense that the drag and drive currents have the same direction, which leads to a 
compensating voltage opposite to the drive current. In the absence of a magnetic field, 
drag phenomena in 2D electron systems are well understood [2j. In that case the drag 
signal can be expressed in terms of the density response functions of each single layer, as 
long as the interlayer coupling is weak Q]. 

Drag experiments have also been performed in the presence of strong magnetic fields, 
where the formation of Landau levels plays a crucial role. Pronounced minima in the 
drag signal were observed at odd integer filling factors for magnetic fields far below the 



regime where spin splitting of Landau levels is seen in the intralayer resistivity 



as 



Completely unexpected was the discovery of a negative drag signal in case of a suitably 
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chosen density mismatch between the two layers la 13. Previous theoretical descriptions 

nfi 

|2J, la] would predict positive drag at any filling factor. In more recent work a mechanism 



for a sign change of drag in strong magnetic fields was found [jj, yielding however negative 
drag for equal densities in the two layers, unlike the experimental situation. 

New hints and constraints for a theory result from a recent detailed experimental 
analysis of the temperature dependence of the drag resistivity by Muraki et al. Q|. At 
high temperatures, ksT ^> huj c , where uj c is the cyclotron frequency, the drag resistivity 
Pd follows the conventional T 2 behavior. It decreases with decreasing T down to scales 
well below huj c , but then, in the data shown for v = 2n + 1 with n — 2, 3, 4, 5, rises again to 
form a pronounced peak at temperatures one order of magnitude below hu c . At still lower 
T the drag resistivity decreases rapidly, consistent with exponentially activated behavior, 
Pr, oc e~ A ^ kBT \ The //-dependence of the activation energy A oscillates with minima 
near half-integer filling factors. The longitudinal intralayer resistivity p xx also exhibits 
activated behavior at the lowest temperature, with an activation energy matching with 
A for those v where both gaps could be extracted from the data. For a density mismatch 
bv = 1 between the two layers a low temperature peak with the same shape as for equal 
densities, but with opposite (that is negative) sign was observed. 

Muraki et al. 11Q] interpreted their data in terms of particle-hole asymmetry and dis- 
order induced localization. Landau levels are broadened by disorder, leading thus to a 
band of energies for each level. Anomalous drag is observed for parameters where the 
Landau level broadening is much smaller than hu c . The longitudinal resistivity p xx of 
each layer as well as the drag resistivity po is due to electrons in extended states of the 
highest occupied Landau level. If the states with energies near the Fermi level ep are 
localized, charge carriers in extended states can be created only by thermal activation 
or scattering across the energy gap between ep and the nearest extended state energy 
level. This explains the observation of activated behavior of p xx and pn at low T. Muraki 
et al. also realized that odd integer filling factors lead to a particle-hole symmetric 
occupation of the partially filled Landau level. Since par ticle-hole asymmetry is known to 
be crucial for drag in the absence of a magnetic field [ljj , it is indeed natural to relate the 
minima in pp observed at odd integer v to this symmetry. Judging from the experiments 
a necessary requirement for negative drag seems to be that the partially occupied Landau 
level in one layer be more, the other less than half-filled. This was noticed already by 
Fmg etal. a.d led tot o speculated elects in a more than half-filled Landau 
level should acquire a hole-like dispersion relation due to disorder. 
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A comprehensive transport theory for the intralayer and drag resistivity in a strong 
magnetic field, which captures localization and extended states, is not available. In this 
work we present a simple semiclassical picture which explains why Landau quantized 
electrons moving in a smooth disorder potential behave effectively like band electrons 
with an electron-like dispersion in the lower and a hole-like dispersion in the upper half of 
a disorder broadened Landau level, and how this leads to the observed anomalous drag. 

Let us first consider the relevant length and energy scales in the double-layer system. 
The distance between the layers varies from 30 to 120 nm. The disorder is due to remote 
donors, which leads to a smooth disorder potential with a correlation length £ of order 
50 — 100 nm. The Landau level broadening, which is related to the amplitude of the 
disorder potential, has been estimated from low-field Shubnikov-de-Haas oscillations to 
lie in the range 0.5-2 K 0, Q]. The magnetic fields at which drag minima at odd integer 
filling and/or negative drag at mismatched densities are observed vary over a relatively 
wide range from 0.1 T to IT for the cleanest samples, and up to 5 T for samples with a 
slightly reduced mobility. This corresponds to magnetic lengths Ib between 15 and 80 nm. 
Hence, £ and Ib are generally of the same order of magnitude, which makes a quantitative 
theoretical analysis rather difficult. Since the anomalous drag phenomena are observed 
over a wide range of fields, one may hope that qualitative insight can be gained also by 
analysing limiting cases. For ^ <C ^ a treatment of disorder within the self-consistent 
3orn approximation is possible This route was taken most recently by Gornyi et al. 
~\\ , who obtained negative contributions to the drag resistivity at mismatched densities in 
agreement with experiment. Localization is not captured by the Born approximation, and 
consequently the resistivities obey power-laws rather than activated behavior for T — > 0. 
For stronger magnetic fields, on the other hand, localization will become important, and 
a semiclassical approximation is a better starting point, which is the route we take here. 
Applying a criterion derived by Fogler et al. we estimate that (classical) localization 
may set in already at 0.1 T. 

We now discuss the semiclassical picture of electron states and drag on a qualitative 
level, ignoring spin splitting for simplicity. If the disorder potential varies smoothly, 
electron states lie essentially on contours of equal energy. These contours form closed 
loops, corresponding to localized states, except at a single energy e in the center of each 
Landau level, for which there is a percolating contour through the whole system jlfil |. 
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If the Landau level broadening induced by the disorder potential and also k B T is much 
smaller than tkv c , as is the case in the anomalous drag regime, all Landau levels except 
one are either fully occupied in the bulk of the whole sample or completely empty. At 
zero temperature and for ep < eo the sample then consists of islands where the highest 
(partially) occupied Landau level is locally full, while it is locally empty in the rest of the 
sample; for tp > eo the empty regions form islands. In reality, the percolating contour at 
the center of the Landau level is broadened to a percolation region for two reasons. First, 
electrons near the percolating contour can hop across saddlepoints from one closed loop to 
another, by using for a moment some of their cyclotron energy [l^ . Second, electrons near 
the center of the Landau level can screen the disorder potential, such that the percolating 
equipotential line at eo broadens to form equipotential terraces Hence, there is a 

region around the percolating path, where states are extended. 

The percolating region forms a two-dimensional random network consisting of links and 
crossing points For simplicity we assume that the links are straight lines. The region 
near the link (except near the end points) can be parametrized by a cartesian coordinate 
system with a variable x following the equipotential lines parallel to the link and y for the 
transverse direction. The disorder potential varies essentially only in transverse direction, 
and can thus be represented by a function U(y). In the Landau gauge A = (—By, 0) the 
Hamiltonian for electrons on the link is then translation invariant in x-direction and the 
Landau states (h = 1 from now on) 

iM*. y) = C n e -(y-^f/VD Hn [(y - l 2 B k)/l B ] e** (1) 

are accurate solutions of the stationary Schrodinger equation. Here H n is the n-th Hermite 
polynomial and C n a normalization constant. The corresponding energy is simply 

e„ fc = (n+ \)u c + U{l B k) . (2) 

In the following we drop the Landau level index n, because only the highest occupied level 
is relevant. The quantum number k is the momentum associated with the translation 
invariance in x-direction. It is proportional to a transverse shift yo = l B k of the wave 
function. The potential U lifts the degeneracy of the Landau levels and makes energies 
depend on the momentum along the link. The group velocity 

v k =^ = l 2 B U'(y ) (3) 
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corresponds to the classical drift velocity of an electron in crossed electric and magnetic 
fields. In our simplified straight link approximation U'(y ) does not depend on x. In 
general it will depend slowly on the longitudinal coordinate, but near percolating paths 
away from crossing points it has a fixed sign, and hence the group velocity a fixed direction, 
that is the motion on the links is chiral. 

The drag between two parallel layers is dominated by electrons in the percolating 
region, corresponding to states near the center of the Landau level, since electrons in 
deeper localized states are not easily dragged along. In the network picture, macroscopic 
drag arises as a sum of contributions from interlayer scattering processes between electrons 
in different links. If no current is imposed in the drive layer, both layers are in thermal 
equilibrium and the currents on the various links cancel each other on average. Now 
assume that a small finite current is switched on in the drive layer such that the electrons 
move predominantly in the direction of the positive x-axis (the electric current moving in 
the opposite direction). This means that the current on links oriented in the positive x- 
direction is typically larger than the current on links oriented in the negative x-direction. 
Interlayer scattering processes lead to momentum transfers between electrons in the drive 
and drag layer. The preferred direction of momentum transfers is such that the scattering 
processes tend to reduce the current in the drive layer, that is the interlayer interaction 
leads to friction. Now the crucial point is that electrons moving in the disorder potential 
are not necessarily accelerated by gaining momentum in the direction of their motion, or 
slowed down by loosing momentum (as for free electrons). The fastest electrons are those 
near the center of the Landau level: they have the highest group velocity on the links of the 
percolation network and they get most easily across the saddlepoints. Electrons in states 
below the Landau level center are thus accelerated by gaining some extra momentum 
in the direction of their motion, but electrons with an energy above e are pushed to 
still higher energy by adding momentum and are thus slowed down. To understand how 
negative drag arises, consider the situation with the highest occupied Landau level of the 
drive layer less than, and the one in the drag layer more than half filled. In that case the 
electrons in the drag layer receive momentum transfers with a predominantly positive x- 
component. As a consequence, electrons near the Fermi level of the drag layer are mostly 
slowed down if they move on links toward positive x-direction, and accelerated if they 
move on links in the negative x-direction, which is the opposite of what free electrons 
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would do. If no net current is allowed to flow in the drag layer, an electric field is built up 
in negative x-direction which compensates the effective force generated by the scattering 
processes with the drive layer electrons. The drag signal is thus negative. 

To substantiate the above qualitative picture, let us analyze the drag between two links 
using semiclassical transport theory. We first discuss parallel links in some detail, and 
then briefly the more general case of a finite angle between the links. For a fixed Landau 
level index, the electronic states in each link are fully labelled by their momentum k. The 
(non-equilibrium) occupation of states in the links is described by a distribution function 
fa{k), where a = 1, 2 labels the two layers. We consider the experimental standard setup 
where a small current flows through the "drive-layer" (a — 2), generating a compensating 
drag voltage in the "drag-layer" (a = 1). No net current is allowed to flow in the drag- 
layer. Under these conditions, the linear response of the drag layer to the drive current 
is determined by the linearized Boltzmann equation 

~dh 



dt 



(4) 



coll 



with k\ = —eEi, where E\ is the electric field leading to the drag voltage. The interlayer 
collision term is given by 

12 



dt 



dfl] =-j^j^ wll krKK iMh) + Mki) - MK) - MK)} 

(* 2 )[1 - - S(e kl + e k2 - e k[ - e k >) , (5) 



coll 



where W^^.^ k , is the rate for a single interlayer scattering event k x — > k[, k 2 — > k' 2 and 
the deviation (from equilibrium) functions ip a are defined as usually by / — f° = — 
/°) ip. There is also an intralayer scattering term with a similar structure. The distribution 
function in the drive layer is obtained from the single layer Boltzmann equation as ip2(k) = 
— e J F2T 3i v k , where j 2 is the drive current and vp 2 the equilibrium Fermi velocity. For 
weak interlayer coupling, the deviation function tpi can be neglected in the interlayer 
collision term, because it would yield a contribution of order (W 12 ) 2 . To determine the 
relation between the drive current j 2 and the drag field Ei, we multiply the Boltzmann 
equation by v kl and integrate over k\. The left hand side yields — e J * fi'(e kl ) v% Ei 
which tends to ^ vfi Ei at low temperatures, where vfi is the velocity at the Fermi level 
of the drive layer. Using the antisymmetry of the integrand of Eq. (5) under exchange of 
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k\ and k[, the right hand side can be written as 
7T/2 f dh f dk 2 f dq TT/12 



e^ 2 T ./ 2tt 7 2tt J 2n klM]kl+q ' k "- q 
,\f?M - me kl+q )\ [#(e fe ) - J 2 °(e fc2 _ ? )] 



5(efei + e fe2 - e fcl+g - e& 2 _ g ) ■ (6) 



4sinh 2 [(e fcl+(? -e fcl )/2T] 
The integrated intralayer scattering contributions cancel due to the condition of vanish- 
ing drag current. Note that the above integral shares several features with the general 
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expression for the drag response function, as obtained from the Kubo formula 
In particular it is symmetric in the layer indices, a property that depends crucially on 
the correct form of ip 2 {k)- For a simplified discussion of the most important points we 
assume that the interlayer scattering rate W 12 depends only on momentum transfers q 
and energy transfers uo = e^+q — e^, and that momentum transfers are so small that one 
can approximate Vk+ q — Vk by q -Jr-. The drag resistivity pn = —Ei/j 2 can then be written 
as 

2vre 2 v F1 v F2 m l m 2 J J-oo sxnh z (uj/2T) 

(7) 

where Xa{q, is the dynamical density correlation function in layer a, and the effective 
masses are given by the curvature of the dispersion relations at the Fermi level 

1 dVka 



m a dk kpo 

The integral in Eq. (7) is always positive. The sign of po is thus given by the sign of the 
effective masses, that is by the curvature of the dispersion at the Fermi level. Negative 
drag is obtained when the dispersion in one layer is electron-like, and hole-like in the 
other. The drag vanishes if the Fermi level in one of the layers is at an inflection point 
of the dispersion. For a quadratic dispersion = J^- one has v F = kp/m and Eq. 
(7) reduces to a one-dimensional version of the well-known semiclassical result for drag 
between free electrons in two dimensions |2J. Returning to Eq. (6), it is not hard to 
generalize the above results on the sign of po allowing for larger momentum transfers q 
and general momentum dependences of W 12 . Note that in our case of chiral electrons no 
backscattering is possible, unlike the situation in quantum wires [20]. 

For parallel links, energy and momentum conservation restrict the allowed scattering 
processes very strongly. At low temperatures, this leads to an exponential suppression 
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of the drag between parallel links. This has nothing to do with the exponential suppres- 
sion of drag observed in the experiments, since the links are generically not parallel. For 
non-parallel links the sum of momenta on the two links is no longer conserved in the scat- 
tering process. Hence scattering processes are suppressed much less at low temperatures. 
Computing the drag between non-parallel links from the linearized Boltzmann equation 
(a straightforward generalization of the above steps for parallel links) yields a quadratic 
temperature dependence at low T. The momentum transfers in the drag and drive links 
are however still correlated for non-parallel links, especially when the angle between the 
links is not very large, and the relative sign of the curvature of the dispersion in drive 
and drag layer, respectively, determines the sign of the drag. The average curvature van- 
ishes for states in the center of the Landau level, while it is positive for energies below 
and negative for energies above 6q. We thus understand the observation of negative drag 
when the Landau level in one layer is less than half-filled, and more than half-filled in the 
other. 

Spin can be easily included in the above picture. Since the interlayer interaction is 
spin independent, one simply has to sum over the two spin species (up and down) in 
both drive and drag layer, taking the (exchange enhanced) Zeeman spin splitting of the 
Landau levels into account. If the Fermi level of one layer lies between the centers of the 
highest occupied Landau levels for up and down spins, respectively, positive and negative 
contributions to the drag partially cancel each other. The cancellation is complete due to 
particle- hole symmetry in the case of odd integer filling, as observed in experiment. 

Within our semiclassical picture anomalous drag, especially negative drag, is sup- 
pressed at temperatures above the Landau level width, because then electron- and hole- 
like states within the highest occupied level are almost equally populated. This agrees 
with the results from the Born approximation Q| , and also with experiments. For the low 
temperature asymptotics of the drag, the semiclassical theory yields two different types 
of behavior, depending on the filling. If the Fermi level does not hit any extended states 
(for either spin species), the drag should vanish exponentially for T — > 0, since thermal 
activation or scattering of electrons into extended states is then suppressed by an energy 
gap. By contrast, for a Fermi level within the extended states band (for at least one spin 
species) the gap vanishes and the drag obeys generally quadratic low temperature behav- 
ior, as obtained for the drag between non-parallel links. Within the Born approximation 
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no localization occurs and the drag resistance always vanishes quadratically in the low 
temperature limit Q| . In high mobility samples localization is negligible at low magnetic 
fields, while an increasing amount of states gets localized at higher fields 

In summary, we have presented a semiclassical theory for electron drag between two 
parallel two-dimensional electron systems in a strong magnetic field, which provides a 



transparent picture of the most sa 
observed in recent experiments 



ient qualitative features of anomalous drag phenomena 
Lo| . Localization plays a role in explaining activated 
low temperature behavior, but is not crucial for anomalous (especially negative) drag per 
se. A quantitative theory of drag which covers the whole range from low magnetic fields, 

nn 

where the Born approximation is valid 13, 14]. to high fields, where localization becomes 
important, remains an important challenge for work to be done in the future. 
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